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>i 

' Gamma-ray bursts (GRBs) are the most brilliant objects in the Universe but efforts 

^ to estimate the total energy released in the explosion — a crucial physical quantity — 

2 have been stymied by their unknown geometry: spheres or cones. We report on a 
O 

comprehensive analysis of GRB afterglows and derive their conical opening angles. We 

^p^nd that the gamma-ray energy release, corrected for geometry, is narrowly clustered 
O 

^ around 5 x 10^° erg. We draw three conclusions. First, the central engines of GRBs 
release energies that are comparable to ordinary supernovae, suggesting a connection. 

> 

Second, the wide variation in fluence and luminosity of GRBs is due entirely to a 
^ distribution of opening angles. Third, only a small fraction of GRBs are visible to a 
given observer and the true GRB rate is at least a factor of 500 times larger than the 
observed rate. 

This paper has been submitted to Nature. You are free to use the results here for the purpose of 
your research. In accordance with the editorial policy of Nature, we request that you not discuss this 
result in the popular press. If you have any question or need clarifications please contact Dale Frail, 
dfrail@nrao.edu or Shri Kulkarni, srk@astro.caltech.edu. 
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Observations of GRBs are well described by the fireball model/ in which an explosive flow of 
relativistic matter (ejecta) is released from a central source. The collision of fast-moving ejecta 
with slower moving ejecta result in bursts of gamma rays. Shortly thereafter, the ejecta starts 
shocking and sweeping up significant amounts of circumburst matter. The shocked gas, hereafter 
called the blast wave, powers long-lived and broad-band (X-ray, optical and radio) emission - the 
so-called afterglow emission. 

The afterglow emission appears to be primarily synchrotron radiation. As with other astro- 
physical shocks, the shocked electrons are accelerated to relativistic energies, forming a power-law 
distribution, dN/dE^ oc E~'P] here, is the energy of the electron and p, the index of the power 
law. In the presence of magnetic fields, the electrons radiate synchrotron emission with a flux 
/(t, z/) oc t°z/'^, where the spectral index (/?) and the temporal index (a) are related to p and the 
dynamics of the blast wave.^ Broad-band observations have repeatedly confirmed the expectations 
of this simple picture. ^"^ 

The two outstanding issues in this field are (a) determining the progenitors of GRBs and (b) un- 
derstanding the physics of the central engine. The focus of this Article is the latter topic, specifically 
the energetics of these mysterious sources. 

Observationally it is known^ that the fluence (defined as the received energy per unit area) of 
the broad-band afterglow phase is always (and usually much) smaller than that of the gamma-ray 
burst fiuence. This then motivates the use of the isotropic equivalent gamma-ray energy, -Eiso(7) = 
47rF^(i|(l + z)~^ as a surrogate for the energy released by the central engine. Here, F-y is the fiuence 
of the burst; z is the redshift; and is the luminosity distance. 

However, £^iso(7) could grossly overstate the true gamma-ray energy release {E^) if the explosion 
is not spherical. Indeed, jets are present in almost all accretion-driven phenomena, e.g., young stellar 
objects, neutron star binaries, and quasars. Likewise, there is excellent observational evidence for 
GRB fireballs with conical geometry. Henceforth, following standard usage, we will interchangeably 
use the term "jet" for conical blast waves. 

As a result of relativistic beaming ( "abberation" ) , an observer can see only a hmited portion of 
the blast wave with angular size ~ F"^, where F is the bulk Lorentz factor of the blast wave. This 
relativistic beaming implies that there is no observable distinction between a spherically expanding 
blast wave and a conical blast wave (whose opening angle we denote by 6j) until the blast wave 
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has slowed down to F < 9J^. However, gamma-ray emission is expected to occur when F is large, 
F ^ 100, and thus unless the opening angles are very small, 9j < 0.01, a conical GRB will not hght 
up the full celestial sphere, but only the so-called beaming fraction = (1 — cos9j); we note that 
for 9j 9]/2. 

There are two consequences of relativistic beaming: (1) the true GRB rate is f^^ times larger 
than the observed GRB rate and (2) the true gamma-ray energy released is smaller than Eigoi'j) 
by the same factor, i.e., = ft, x Eisoi'j)- 

In contrast to the situation during the gamma-ray burst phase, during the afterglow phase F 
ranges from ~ 10 hours after a burst to trans-relativistic (F ~ 1) and even non-relativistic values 
days to months after the burst. Thus at some point during the lifetime of the afterglow F falls below 
9~^ which has a clear observational signature. Thus mult i- wavelength (X-ray, optical and radio) 
afterglow observations offer us an elegant way to measure 9j. Here, we report an extensive analysis 
of 9j values for all well-studied GRBs. Prom these values of 9j we are able to infer the true 7-ray 
energy release of the central engines of GRBs. 

The organization of the paper is as follows. First wc summarize the physics of conical afterglows 
followed by extraction of 9j values from afterglow observations. We then find the surprising result 
that E^ is tightly clustered around 5 x lO^'^ erg. We end with a discussion of the ramifications of 
this result. 

Conical Afterglows 

The temporal evolution of the afterglow emission is directly related to the dynamics of the blast 
wave which in turn is influenced both by the circumburst medium and the geometry of the explosion. 
The afterglow emission from a conical blast wave^"^^ differs from that of a spherical blast wave in 
two distinct ways. First, the observer will start noticing a deficit of emitting material when F < 9J^. 
The magnitude of this deficit, relative to that of a spherical fireball, is proportional to the ratio of 
the area of the emitting surface for a conical blast wave (oc 9j) to that of a spherically emitting 
blast wave (oc F~^). This deficit results in the afterglow emission declining more rapidly, relative to 
a spherical case, or a "break" in the power law decay, Aa = 3/4. 

The second effect that also becomes important when F ^ 9~^ is the spreading of the jet in the 
lateral dimension. The ejecta now encounter more surrounding matter and decelerate faster 
than in the spherical case. This results in an overall steepening (Act — 1.2-1.5) of the afterglow 
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emission. In the case of a laterally spreading jet, the lightcurve evolves with a = p. We note p ~ 
2.2-2.4 appears to fit all well-studied afterglows (e.g., ref. 11, 12). This value is also favored by shock 
acceleration models. 

The first claim of a jet was made for the radio afterglow of GRB 970508, which showed deviations 
from the predictions of a simple spherical adiabatic model. ^ However, it was the spectacular isotropic 
energy release^^ of GRB 990123 - approaching the rest mass of a neutron star - which emphasized 
the possible importance of jets in GRBs. A case for a jet in the afterglow of this burst was made on 
the basis of a sharp break (Aa > 0.7)^^ in the optical afterglow and upper limits in the radio. The 
clearest evidence for a jet is a sharp break over a broad range of frequencies and such a signature 
was seen in the hghtcurves of GRB 990510 at optical^^'^^ and radio^^ wavelengths and was found 
to be consistent with the X-ray^^ light curve. Furthermore, the detection of polarization^^'^^ from 
this event gave further credence to the jet hypothesis: the non-spherical geometry leads to polarized 
signal, from which the geometry of the jet can be inferred. ^^'^^ 

More recently, the identification of jets has shifted from single frequency measurements to global 
model fitting of joint optical, radio and X-ray datascts (e.g., ref. 5, 24). This approach has the 
advantage that by simultaneously fitting all the data, the final outcome is less sensitive to deviations 
in small subsets of the data. In addition, since the character of the achromatic break is different 
above and below the peak of the synchrotron spectrum, broad-band measurements give more 
robust determinations of the jet parameters. This approach was crucial in distinguishing the jet 
break for GRB 00030 IC^ whose decaying lightcurves exhibited unusual variability,^^ now attributed 
to microlensing.^^ 

Determination of the Jet Opening Angles 

We use the formulation of Sari, Piran & Halpern^^ to convert the measured jet break times tj to 
opening angles of the conical blast wave: 

^ Viday^ V 2 ^ Vio53erg/ vo.2^ VQ.lcm-s/ ' ^ ' 

where 77-y is the efficiency of the fireball in converting the energy in the ejecta into 7 rays, and n is 
the mean circumburst density. In Table 1 we present a complete sample of all GRBs with known 
redshifts as of December 2000. The determinations of tj are of varying quality. The best events are 
those for which it is possible to globally model the broad-band data within the physical framework of 
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the relativistic jet model (e.g., GRB 000301C, GRB 990510). For some bursts tj is inferred from only 
one band (e.g., GRB 990705) and in some cases with additional constraints from radio observations 
(e.g., GRB 990123). Finally, there are some events with only upper (e.g., GRB 991208) or lower 
limits on tj (e.g., GRB 971214), for which only upper or lower limits of 9j can be placed, respectively. 

We obtain a range in 9j corresponding to the wide range in tj values in Table 1 (from ^ 1 d to 
30 d). The derived jet angles vary from 3° to more than 25° with a strong concentration near 4° 
(Figure 1). It is reasonable to ask whether the observed distribution in Figure 1 suffers from selection 
effects. To begin we note that out of the 21 known optical afterglows, the light curves of only two 
GRBs - GRB 980326 (ref. 27) and GRB 980519 (ref. 28) - show rapid dechne implying tj ^ 1 d. 
Likewise, out of a sample of 10 bright X-ray afterglows observed with the BeppoSAX sateUite there 
is no evidence for a significant break within 8 to 48 hours after a burst, suggesting that tj ^ Id 
for these events. If we increase the sample to include the 28 GRBs detected by BeppoSAX for 
which follow-up searches (typically 8-12 hr after the burst) were made for an X-ray afterglow we 
find only one unambiguous case where no afterglow was detected (GRB 990217; ref. 30). There are 
a further six cases where a hitherto uncataloged X-ray source was detected in the GRB error circle. 
In every case the X-ray source is a plausible afterglow but lacking multi- wavelength confirmation, 
the afterglow identification remains uncertain, e.g., GRB 970111, ref. 31. From these statistics we 
conclude that steep decays, tj ^ 1 d, and therefore very narrow opening angles, 9j < 3°, are required 
for less than ten percent of the BeppoSAX GRB sample. 

There is another method to infer the existence of a population of GRBs with extremely narrow 
opening angles. The beaming fraction during the afterglow phase is max(^|/2, r~^/2). Thus, while 
narrow-angle GRBs will be rare, their X-ray, optical, and radio afterglows which are emitted at 
increasingly smaller F are accordingly less rare.^ However the current limits^^'^^ of these "burst- less" 
afterglows do not place further significant constraints on 9j. 

GRBs with large opening angles do not suffer from severe beaming but it is not easy to measure 
tj for such bursts. For large tj the afterglow emission is weak and (at optical wavelengths) the 
host galaxy starts dominating. ^^'^^ Thus optical observations and X-ray are unhkely to yield tj. 
Fortunately, radio observations can and do play a crucial role, due to the long lifetime of the afterglow 
in this regime. This was the case for four out of five wide-angle jets identified in Table 1. One notable 
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example is GRB 970508 where a jet modeP^ of the radio data was found to be consistent with an 
analysis of the optical hght curves. ^° 

The Luminosities and Energies of GRB Central Engines 

In Figure 2 we plot the measured fiuence versus the inferred inverse beaming factor. There 
appears to be a correlation in the sense that the bursts with the largest fiuence have the narrowest 
opening angles. This trend was noted earher^^ albeit based on a few afterglows. The correlation is 
improved when the fluences are all scaled to the same redshift (unity), which effectively renders it 
to a correlation between Eisoi'j) and fj^^. The physical meaning of this trend is better appreciated 
from Figure 3 where we find that E^, the true energy released in gamma rays, is clustered around 
5 X 10^° erg, with a la multiplicative factor of only two. 

Figures 2 and 3 suggest the following simple scenario: the central engines of GRBs produce 
approximately a similar amount of energy, and a significant part, about 10^^ erg, escapes as gamma- 
rays (Figure 3). However, for reasons not presently understood, there exists a wide range of jet 
opening angles. If so, GRBs with the narrowest opening angles would be brighter and consequently 
produce the correlation seen in Figure 2. 

The narrowness of the distribution is surprising and has several immediate imphcations. While 
it is not unreasonable to expect that the central engines produce a similar amount of energy, Eq, in 
each explosion, there is little reason to expect that they will produce similar gamma-ray outputs. 
Since the true total energy Eq = E^/rj^ oc n^/^jy"^/^ (this follows from Equation 1), the narrowness 
in the distribution of E^ places restrictions on the dispersion of n and r^^. 

If ?7-y is high (close to unity) then a small dispersion in r]^ is naturally assured. Indeed, a number 
of recent papers^^"^^ have argued that internal shocks under certain conditions are very efficient at 
producing gamma rays {rjj ^ 0.2). Furthermore, Guetta et al.^^ argue that the very conditions that 
are needed to make internal shocks efficient (a large dispersion in the distribution of the ejecta's 
Lorentz factors) also produce the characteristic clustering of spectral break energies of GRBs in the 
range 0.1-1 MeV. 

Given a distribution in E^ with a full width of a factor of four (see Figure 3) , the dispersion in n 
{(X E^) has to be less than two orders of magnitude. At first this may seem to be a weak constraint on 
the possible progenitors. However, the progenitors discussed to date lie either in intergalactic space 
or the halos of galaxies (ns-ns coalescence, n ~ 10~^ cm~^ and n ~ 10~^ cm~^) or in a typical disk 
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interstellar medium (ISM) or dense ISM (coUapsar; n ~ 1 cm~^ and n ~ 10^ cm~^, respectively). 
Therefore our results limit the diversity of GRB environments, and specifically requires that the 
long-duration class of GRB events happen in only one of these environments. Furthermore, we note 
that winds of massive stars would produce a density of a few atoms cm~^ for tj 1 d. In scenarios 
where there are two types of GRBs,^° the ones that do not go off in stellar- wind-stratified media must 
reside within the disk of their host galaxy rather than in galaxy halos or the intergalactic medium. 
Indeed, the distribution of GRBs within their host galaxies is consistent with a disk population.''^ 
Likewise, broad-band modeling of GRB afterglows^'^^'^^ give estimates of gas densities consistent 
with disks, justifying our normalization of n in Equation 1. We conclude that the progenitors of 
long duration GRBs likely come from one type of progenitor. 

Finally, the narrowness of the Ej distribution requires that the brightness of the 7-ray beam be 
roughly uniform from the center to the edge. This is contrary to models^^ in which large intensity 
variations within the conical blast wave are invoked in order to explain the wide dispersion of peak 
luminosities. We find that most of the dispersion in the luminosity is due to the diversity in opening 
angles. 

The mean value of £^-y is 5 x 10^° erg (Figure 3). If we accept the conclusions of Guetta et al.^^ 
(see above), then 77^ ~ 0.2 and we then derive £'0 ~ 3 x 10^^ erg. Of course, Eq is sensitive, in 
addition to the adopted value of r)^, to the overall scaling, i.e., the numerical coefficient of Equation 
1. For example, the estimate of Rhoads,^° based on a different assumption for the sideways expansion 
speed, has a coefficient smaller by a factor of six than our Equation 1. 

Fortunately, GRB 970508 allows us to directly determine the energy scale. The radio afterglow 
of this GRB lasted long enough (400 d) that the blast wave was non-relativistic, thereby allowing 
determination of the total energy^^ independent of relativistic beaming. Table 1 shows that this 
burst has one of the lowest energies, although it is only la away from the mean (if the energy 
distribution is assumed to be log normal). The agreement between these two entirely different 
approaches is remarkably good and gives some support to our choice of the numerical coefficient and 
normalization of Equation 1. 

Freedman & Waxman^^ and Kumar"'^ have suggested an elegant way to estimate the energy 
in the afterglow phase based on X-ray observations. This method yields the e^ea where is the 
energy of the blast wave per steradian and eg is the fraction of energy in the shocked electrons. This 
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estimate is independent of the ambient density. If is high and relatively constant (analogous to the 
situation with 77^) then Eq can be estimated provided /{, is known. It is of interest to note that the 
ratio Ea/Eisoi'j) is nearly constant/^ suggesting that likely both eg and rj^ are narrowly distributed. 

Applying our determinations of fh to the sample of ref. 12 (a total of six common GRBs) we 
obtain Ea — (2.7 ± 1.4) x 10^*^ erg. Within the limitations of the small sample, the distribution 
appears to be clustered and the results are in agreement with our findings (Figure 3). The principal 
advantage of our method is that the events are always identified in the 7-ray band, whereas X-ray 
observations are available for only a minority of cases. Furthermore, the X-ray afterglow technique 
ignores the effects of inverse Compton scattering effects, and is therefore sensitive to the poorly 
known strengths of magnetic fields in strong shocks. 

GRBs and SNe 

Above we find that the mean total energy of GRBs is £^0 ~ 3 x 10^^ erg. This energy is only 
slightly larger than the typical 10^^ erg of electromagnetic and kinetic energy yield of ordinary 
supernovae (la, Ibc, II). This reduced energy budget raises the possibility that GRBs are the result 
of the formation of neutron stars,^^ albeit with special properties, and does not necessarily require 
black holes. The mystery about GRBs is no longer in understanding their supposedly extraordinary 
energy budget but in explaining why the ejecta of GRBs have such a high Lorentz factor. 

We note however that there are at least two possible exceptions to the tight clustering of jet 
energy. (1) If SN 1998bw is associated with a GRB^^'^^ then Ei^{-f) ~ 7 x 10^^ erg. However, 
Kulkarni et al.^^ have argued that the extraordinarily bright radio emission from this SN requires 

]^q50 gj,g ^j^g explosion energy to be in the form of mildly relativistic ejecta (F ~ few). (2) Bloom 
et al."^^ identify the late time red bump in the rapidly decaying event GRB 980326 {tj < 0.55 d; ref. 
27) with an underlying SN. If so, the inferred redshift 2; ~ 1 and E^ < 7 x 10^^ erg. Unfortunately, 
the radio observations are not sensitive enough to place meaningful constraints on the amount of 
energy in mildly relativistic ejecta. In both cases, the true energy release could be closer to Eq but 
this energy could be primarily in mildly relativistic ejecta. Careful observations (especially X-ray 
and radio) of SNe may uncover significant numbers of such "failed" GRBs. 

Beaming Fraction and the GRB Rate 

Since conical fireballs are visible to only a fraction, /b, of observers, the true GRB rate, 
— {fb^)Robs: where Rohs is the observed GRB rate and (f^^) is the harmonic mean of the 
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beaming fractions. We find (f^^) ~ 500 (see caption to Figure 1). The formal uncertainty in this 
estimate is only 16% but systematic uncertainties related our choice of the numerical coefficient and 
normalization of Equation 1. make this estimate accurate to a factor of two. 

Estimates^°"^^ of the local observed rate of GRBs give values of Robs{z — 0) ranging from 0.2 
to 0.7 Gpc"^ yr~^. The rate is uncertain because it is not known how the GRB rate evolves 
with redshift. We adopt a value -Robs(-2^ = 0) = 0.5 Gpc~^ yr~^ as in Ref. 51. The true rate is 
Rt{z = 0) ~ 250 Gpc~^ yr~^, which should be compared with the estimated rate^^ of neutron star 
coalescence, Rc{z = 0) ~ 80 Gpc~^ yr~^ and the estimated rate^^ of type Ibc SN, i?ibc ~ 6 x 10^ 
Gpc~^ yr~^. Clearly, the coUapsar scenario is capable of easily supplying a sufficient number of 
progenitors (including failed GRBs). Within the uncertainties of the estimates, the coalescence 
scenario is also (barely) capable of providing sufficient progenitors. 

Assumptions, Uncertainties and Caveats 

Our derivation of the jet opening angle is based on Equation 1 which makes two implicit assump- 
tions. First, we assume that GRBs explode in a constant density medium and that any sharp break 
in the afterglow (Aa 0.75) is attributed to a combination of the observer viewing beyond the edge 
of the conical jet and sideways expansion. Second, we assume that the conical blast wave maintains 
a fixed opening angle right from the GRB phase until F approaches 9^^. The latter assumption 
imposes strict conditions on the working of the central engine. Specifically, the ejecta have to be 
approximately uniform across the entire opening angle in the gamma-ray phase and the bulk of the 
explosive energy in the afterglow phase must have a single bulk Lorentz factor. 

The origin of the observed break is currently a matter of considerable theoretical debate. ^"^"^^ The 
uncertainty is driven by the as yet unclear hydrodynamics of sideways expansion. Some authors^^'^^ 
argue that transition is very smooth, and is completed in one decade of the break time for a constant 
density circumburst medium but takes two decades for the 71(7 ) oc r"^ circumburst medium. Others-*^^ 
argue that several uncertainties in these calculations make this conclusion premature. 

Furthermore, the analysis of the afterglow lightcurves does not yield 0, the angle between the line 
of sight to the observer and the principal axis of the jet. In general one may expect some dependence 
of tj on (f> and thus distribution should broaden even if £"0, rj^ and n were constants. 

However, the narrowness of the distribution shown in Figure 3 provides empirical support for 
our assumptions. Furthermore, as noted earlier, where high quality observations are available (e.g.. 
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GRB 990510 and GRB 000301C) the breaks are found to be quite sharp and the inferred ambient 
density is ~ 1 cm~^. 

Finally, several other mechanisms have been proposed to produce steep declines in the afterglow 
light curves: (i) Si sudden drop in the external density,^^ (ii) transition from relativistic to non- 
relativistic regime^^ due to expansion in a dense circumburst medium, and (in) a break in the 
power-law distribution of radiating electrons. These models have not been systematically compared 
against well studied afterglows and thus remain at a level of suggestions. Model (ii) can be rejected 
because the expected centimeter wave attenuation due to free-free absorption is not seen. Model 
(Hi) is unable^^ to explain the broad-band achromatic breaks and will therefore fail to account for 
the early time low frequency emission. We conclude that at this stage the simple jet model which 
we have adopted provides a consistent and adequate description of the observations. 

Our understanding of gamma-ray bursts has increased dramatically over the past four years. For 
nearly three decades these objects were considered so enigmatic that expectations of their distance 
ranged from local to cosmological scales. In the BATSE era, prior to the discovery of the after- 
glow phenomenon, the standard assumption was that GRBs possessed fixed peak himinositics. As 
more and more redshifts were obtained, the isotropic equivalent energy record increased, eventually 
reaching the rest mass energy of neutron stars. The standard candle hypothesis was consequently 
abandoned. It is remarkable that with a more detailed understanding of the afterglow we are able 
to infer the energy release in these bursts and find that GRBs are "standard candles" in some sense. 

We have deduced the distribution of the opening angles of GRB jets and empirically uncovered 
a key clue, namely the total energy release and its approximate constancy, but we are still left 
with three significant mysteries. First, we do not know what physical mechanism results in such 
a wide variation in the opening angles of the jets. Second, the similarity of the energy release in 
GRBs and ordinary supernovae is puzzling. This coincidence is all the more remarkable considering 
the diversity of the progenitors and likely differing collapse mechanisms in these various classes of 
explosions. Third, we do not understand why in GRBs, the explosion energy couples only to 10~^ Mq 
of the exploding star and thereby produce ejecta with high Lorentz factor. Fortunately, new missions 
(HETE-2 and SWIFT) with their vastly increased GRB localization rates will provide empirical data 
which may help solve these mysteries. 
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Opening Angle {6., rad) 



Figure 1. The observed distribution of jet opening angles along with a model fit (line). We assume 
that the observed differential distribution of beaming factors can be represented by two power laws: 

Pohsifb) = (fb/fo)"^^ for fb < fo and Pohsifb) = Ub/ kY""^ for h > fo- Since for every observed 
burst there are f^^ that are not observed, the true distribution is ptmeifb) = ft ^Pohsifb)- Fitting to 
the data, we find the following: a is poorly constrained; f3 = — 2.77lo;3o; fog/o = ~2.9llo;og. Thus, 
the true differential probability distribution (under the small angle approximation, fj, (xO^) is given 
by Ptrue(^j) ^7^'^^ with the observed distribution being Po^s ^7^'^^- The distribution Ptmeifb) 
allows us to estimate the true correction factor, {f^^) that has to be apphed to the observed GRB 
rate in order to obtain the true GRB rate. We find (f^^) = fo^P - ~ 520 ± 85. 
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Figure 2. The gamma-ray fluence (in units of erg cm~^) plotted as a function of the inverse 
beaming fraction f^^ (where /(, = ! — cos^j = ^|/2, and 9j is the opening angle of the jet). A 
correlation is apparent in the sense that GRBs that have narrower jet opening angles are brighter 
(high fluence) than those that do not. A linear fit to these data (open squares) gives a relatively 
large rms scatter of a factor of 3.3. The correlation is improved when the fluences are all scaled 
to the same unity redshift (crosses), thereby removing the distance dependence. The rms scatter 
(dashed lines) of these points around a linear fit (solid line) is reduced to only a factor of 2.3. This 
factor is marked by dashed lines around the linear fit. The wide variation in observed fluence, more 
than two orders of magnitude, appears to be mainly due to different beaming angles. 



20 




1049 IQ^O IQ51 lQb2 IQ53 

Energy (erg) 

Figure 3. The distribution of the apparent isotropic 7-ray burst energy of GRBs with known 
redshifts (top) versus the geometry-corrected energy for those GRBs whose afterglows exhibit the 
signature of a non-isotropic outflow (bottom). The mean isotropic equivalent energy {Eisoij)) for 
17 GRBs is 110 X 10^^ erg with a l-cr spreading of a multiphcative factor of 6.2. In estimating the 
mean geometry-corrected energy (E^) we applied the Bayesian inference formalism^" and modified 
to handle datasets containing upper and lower limits. Arrows are plotted for five GRBs to indicate 
upper or lower limits to the geometry-corrected energy. The value of (logE^) is 50.71±0.10 (la) or 
equivalently, the mean geometry-corrected energy (E^) for 15 GRBs is 0.5 x 10^^ erg. The standard 
deviation in log£'-y is O.Sllo oe, or a 1-a spread corresponding to a multiplicative factor of 2.0. 
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Table 1. Jet Break Times and Energetics. The gamma-ray fluences (F^), given in units of 10~^ 

erg cm~^, are from a diverse collection of instruments. The best determinations of energy fluence are from 
the Burst and Transient Experiment (BATSE) on the Compton Gamma-Ray Observatory (CGRO). Most 
of the GRBs (10), prior to the de-orbit of CGRO on 2000 May 26, are BATSE bursts. In these cases 
we used fits to BATSE data which were integrated over the energy range from 20 to 2000 keV. For the 
remainder of the events we used the fluence as determined from the Gamma-Ray Burst Monitor on the 
BeppoSAX satellite (40-700 keV), or fluences (25-1000 keV or 25-100 keV) from the Interplanetary Network 
of satellites {Ulysses, KONUS, and NEAR)7^ The luminosity distance (cIl) is given in units of 10^^ cm. 
It was calculated from the observed redshift (z), and adopting cosmological parameters of Ho=65 km s~^ 
Mpc~^, ^lM=0■3, and Ao=0.7. Other realistic cosmologies were tried but they did not fundamentally change 
our conclusions. The isotropic 7-ray energies (£^iso(7)), given in units of lO^-*^ erg, have been "k-corrected" 
such that all energy estimates are referenced to the same 20-2000 keV co-moving bandpass. Although 
these order-of- unity corrections affect individual determinations of £^180(7), they do not affect our results 
derived from the sample as a whole. The jet break times (tj), given in days, are taken from the literature. 
The notes and the references in the table indicate how tj was determined. The strongest evidence for 
collimated outflows come from GRBs with achromatic breaks in their broad-band light curves (B). In 
most cases such multi-frequency datasets are not available, so there is a second class of events with breaks 
determined primarily from radio (R), optical (O), or X-ray (X) data. We include here a number of events 
for which no break was observed, yielding only lower limits of tj. For some GRBs the steep decline of the 
light curve, indicating a jet geometry, is already fully manifest at the time of the first measurement. In 
these cases (D) we have only an upper limit on tj. The final group of GRBs are those for which tj cannot 
be determined (N) , owing to complications in the light curve such as the presence of a supernova signature 
(i.e., GRB 970828), or the lack of sufficient data. The beaming-corrected gamma-ray energy {Ej), given in 
units of 10^^ erg, was calculated by applying the geometric correction factor to £^iso(7)- 



